Long-range Rydberg-Rydberg interactions in calcium, strontium and ytterbium 
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Long-range dipole-dipole and quadrupole-quadrupole interactions between pairs of Rydberg atoms 
are calculated perturbatively for calcium, strontium and ytterbium within the Coulomb approxi- 
mation. Quantum defects, obtained by fitting existing laser spectroscopic data, are provided for all 
5*, P, D and F series of strontium and for the 3 P2 series of calcium. The results show qualitative 
differences with the alkali metal atoms, including isotropically attractive interactions of the stron- 
tium 5*o states and a greater rarity of Forster resonances. Only two such resonances are identified, 
both in triplet series of strontium. The angular dependence of the long range interaction is briefly 
discussed. 

PACS numbers: 34.20.Cf, 32.80.Ee 
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I. INTRODUCTION 

In the context of laser cooling, there has been a resur- 
gence of interest in the strong interactions between Ry- 
dberg atoms. A key development was the concept of the 
dipole blockade, where the strong interaction between 
Rydberg atoms leads to the excitation of collective states 
with a single Rydberg excitation shared between many 
atoms [l(. Applications of the dipole blockade include 
quantum information [2J], where two-qubit gates have 
already been demonstrated [3|, |j|, co-operative nonlin- 
ear optics Q and the physics of strongly correlated sys- 
tems [a-Q . Other important areas of study in cold Ryd- 
be rg g ases include the formation of long-range molecules 
[9|, [lfj , and the interplay between Rydberg gases and ul- 
tracold plasmas \U\ . 

Central to all of these applications is a detailed un- 
derstanding of the long-range interactions between Ry- 
dberg atoms. As most experiments so far have been 
carried out using the alkali metals (Rb,Cs), theoretical 
work has concentrated on these elements J12| - tl7l ]. Both 
perturbative calculations of t he q uadrupole-quadrupole 
and dipole-dipole interactions |12h14| and detailed non- 
perturbative calculations of the dipole-dipole interactions 
[151 - 1171 have been performed. 

Atoms with two valence electrons offer a new approach, 
as the inner valence electron provides an additional way 
to probe and manipulate Rydberg atoms. Recent experi- 
ments have shown that the inner electron can be used as 
a fast, state-selective probe of the interactions in a cold 
Rydberg gas [18|, [l9| • The polarizability of the additional 
electron also enables tight, magic- wavelength traps for 
Rydberg atoms [20J. Rydberg states of Sr and Yb have 
also been proposed for high-precision measurements of 
the black-body shift in optical frequency standards [21| . 

Preliminary calculations of the interactions for stron- 
tium revealed interesting features that differ from the 
alkali metals, such as the possibility of an isotropic, at- 
tractive interaction potential [2(| ■ The goal of this work 
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is to carry out a systematic study of the interactions be- 
tween Rydberg atoms in the most commonly laser-cooled 
two-electron systems, i.e. Ca, Sr and Yb. In section UU 
we present the background theory for the calculations. 
The C5 and Cq coefficients governing the strength of the 
long-range interactions are expressed in terms of angu- 
lar factors and radial matrix elements. The latter are 
evaluated using the Coulomb approximation, which re- 
quires accurate knowledge of the Rydberg energy levels. 
We therefore present a detailed review of the available 
experimental energy levels. The main results of this pa- 
per are C5 and Cq coefficients for all the Rydberg se- 
ries where experimental energy levels are available. Ta- 
bles of our complete results are provided in the supple- 
mentary data. In section IIIII we present an overview of 
these tables, revealing significant differences between the 
species, and between different spin symmetries (singlet or 
triplet). Several Forster resonances, namely accidental 
near-degeneracies in pair states of the two atoms causing 
a large value of the Cq coefficient, are also identified. 



II. METHODS 

A. Theory 

We consider two interacting divalent atoms, one in the 
state niLiS\J\ and the other in the a state riiLi2S2Ji-, 
separated by a distance R. Specifically, we consider 
states with high enough principal quantum numbers ri\ 
and n% that effects arising from configuration mixing can 
be neglected. Such effects, e.g., singlet-triplet mixing in 
the hsnd X ^D 2 states of strontium [23, can be large in 
the vicinity of perturbers. However, in the species con- 
sidered sharp deviations are only observed close to the 
perturber's energies. When far away from any perturber 
the Rydberg series varies monotonically, thus the effect of 
a perturber can be incorporated into the quantum defect. 
For strontium the perturbers are found much lower in the 
Rydberg series than the states considered here, however 
two of the Rydberg series of calcium and ytterbium are 
considerably perturbed over a wide range of states. A 
single active electron treatment is questionable for the 



calculating the C5 and Cq coefficients affected by these 
perturbations. The effect of these configuration interac- 
tions will be considered in future work. Here, we describe 
each atom by a simple model in which one of the two va- 
lence electrons is in an extended Rydberg orbital while 
the other is in a compact inner orbital. Denoting by r a 
and rfc the position vectors of the valence electrons rela- 
tive to the nucleus, we represent the state of each atom 
by the symmetrized wave function 

Ev^ v^ ^jmj r LM L r iSM s 

Ms-,M L mi,mi m Si ,m Bo 

X [<M r a)Xm s .(a)0o(r6)Xm 3o (>) 

+ (-l) 5 <Mr fc ) XmSi Wo(r a )x ms >)] ■ C 1 ) 

The subscripts i and o refer to the inner and outer or- 
bitals, Xm s (a) and Xm s (b) are spinors describing the spin 
states of electrons a and b, <j>i (r) is the wave function of 
the inner electron, and o (r) the wave function of the 
outer electron. We assume that the inner orbital is of 
s-symmetry; therefore k = m^ = and l = L. 

The Hamiltonian for this system can be written as 



H = H + H s 



int j 



(2) 



where Ho is the Hamiltonian of the pair of atoms at in- 
finite separation and H ln t is the interaction energy be- 
tween the two atoms. In terms of the position vectors of 
the valence electrons and the relative position R of the 
two nuclei 
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with the indexes 1 and 2 identifying the respective atoms, 
and R — |R|. (Atomic units are used throughout this 
section). 

We treat the interaction Hamiltonian H m t perturba- 
tively. For long range Rydberg-Rydberg interactions, the 
matrix elements of H[ nt between the bound eigenstates 
of Hq are dominated by the contribution of the two outer 
electrons. We can therefore reduce the calculation of the 
dispersion coefficients to a problem in which each atom 
has only one active electron and the interaction Hamil- 
tonian reduces to 
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(4) 



where ri and r 2 denote the position vectors of the two 
Rydberg electrons respective to the corresponding nuclei. 
(Although each atom has only one active electron, the 
quantum numbers of the system remain those pertaining 
to the original multi-electron problem.) This formulation 



neglects exchange interactions, which is valid as long as 
R much exceeds the LeRoy radius |23j |. 



i?LR = 2(<r^ + (r 2 2 )^ 



(5) 



A similar single active electron treatment has been previ- 
ously shown to yield accurate Stark maps for two-electron 
Rydberg atoms dHH. 

Expanding H[ nt in multipoles yields [25j 



-"int 
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(47r) 3 (2fei+2fc 2 )! 



y (2fci + l)!(2fc 2 + l)!(2fei + 2fc 2 + 1) 

ki + k 2 fci k 2 

Y v^ V V r kl+k2 ' p 

p=-(kl+k 2 ) Pl =— fcl P2 = -fe2 

x rt 1 r^n ijPl (f 1 )Y fe2 , P2 (f 2 )y fcl+fe2)P (A), (6) 

where R is the unit vector along the internuclear axis. For 
infinite atomic separation, the eigenenergies of the Hamil- 
tonian H coincide with those of Hq, which are sums of 
energies of unperturbed atomic states. These asymptotic 
eigenenergies are thus degenerate in Mi and Af 2 . (There 
is no degeneracy if J\ = ./ 2 = since in this case the 
magnetic quantum numbers M\ and M 2 can take only 
one value.) H[ nt mixes states of different Mj values and 
as a result splits the degenerate asymptotic energy levels 
into a number of sublevels. Each of the latter differs from 
its R — > 00 limit by a sublevel-specific, independent shift 
AE. Treating H[ nt perturbatively and making use of the 
multipolar expansion ([5]) leads to an expression of these 
shifts in the form of a series of inverse powers of R, 



"^1 



N 



(7) 



For the systems we are concerned with, this expansion is 
dominated at large interatomic separations by the term 
in 1/R 5 , when C5 ^ 0, which arises to first order in H' int 
from the quadrupole-quadrupole interaction (k\ = fc 2 = 
2 in the multipolar expansion of H[ nt ), and the term in 
1/R 6 , which arises to second order from the dipole-dipole 
interaction (fci = fc 2 = 1). 

The C5 and Cq coefficients are the eigenvalues of the 
(2Ji + 1)(2J 2 + 1) by (2Ji + 1)(2J 2 + 1) matrices C 5 (-R) 
and Cq(R) formed by the Mj-dependent, /j-dependcnt 
coefficients 

c 5 (M[M! 2 ,M 1 M 2 ;R) = 

D 22 {aM' l M! 2 , aM 1 M 2 ;R)R 22 {n 1 n 1 a, n 2 n 2 a) (8) 

and 

c 6 (M[M^, M 1 M 2 ; R) = - ^ i 

xD 11 (aM[M 2 , ,a"M['M%; R)R 11 (mn 2 a, n"n 2 a") 
xD 11 (a"M['M 2 y ,aM 1 M 2 ;R)R 11 (n'{n' 2 W',n 1 n 2 a). (9) 



In these two equations, a denotes the sextuple of quan- 
tum numbers L1S1J1L2S2J2 and 

Rk 1 k 2 (nin 2 a,n' 1 n' 2 a') = 
Jo 



x / dr 2 Pn 2 L 2 s 2 J 2 (r2)r2 2 Pn'L'S'jd r 2)- 



(10) 



The functions Dn and -D22 are defined in the appendix. 
Moreover, in equation ((9|) the sum runs over all the inter- 
mediate n'{ri 2 a" M" M 2 pair states dipole coupled both 
to the n\ri20iM\M2 state and to the nin2<xM[M 2 state, 
and A denotes the Forster defect (A = E'{ + E 2 ' -E\- 
E 2 ). 

It follows from equation ^ that H[ nt only couples pair 
states of same value of M\ + M2 , and that the ensuing en- 
ergy shifts do not depend on the overall sign of M\ + M2 , 
when the angle 9 between the interatomic axis and the 
axis of quantization of the angular momenta (which we 
take to be the z-axis) is zero. The sublevels the asymp- 
totic energy levels split into are thus characterized by 
|fi|, where fl = M\ + M 2 - ($7 can be recognized as 
the magnetic quantum number associated with J tZl the 
z-component of the total angular momentum operator 
Jt = Ji + J2O For other orientations of the interatomic 
axis, H[ nt also couples pair states differing in il. How- 
ever, the choice of the quantization axis is arbitrary and 
thus the eigenenergies of the system do not depend on R: 
changing the orientation of the interatomic axis changes 
the composition of the eigenstates of the Hamiltonian in 
terms of the unperturbed pair states but does not affect 
the energy shifts A.E. Therefore the dispersion coeffi- 
cients as defined by equation (0 do not depend on R and 
can be obtained by diagonalizing the matrices C§(R) and 
Cg(i?) for any orientation of the interatomic axis. 

Due to the selection rule mentioned at the beginning 
of the last paragraph, these two matrices are block di- 
agonal when this axis is in the z-direction, each block 
being formed by pair states with the same value of ft. 
It follows from the relationships between fl, J t , Ji and 
J 2 that the linear size of each block, i.e., the number of 
values of C5 or C$ associated with each value of il, is 
J\ + J 2 — max( |fi| , I J\ — J 2 1) + 1. In particular, the diago- 
nal blocks with fi = ±(Ji + J 2) contain only one element. 
For these values of O, the dispersion coefficients are thus 
given directly by equations ((8]) and ((9]) as 

C 5 = c 5 (J 1 J 2 ,J 1 J 2 ;R = z) (11) 

C 6 = c 6 (J 1 J 2 ,J 1 J 2 ;R = z). (12) 

($7 = ±( J\ + J2) implies that each atom is in a stretched 
state with M x = ±J X and M 2 = ±J 2 .) 

The calculation thus amounts to diagonalizing the 
Hamiltonians 



H^ = H 



E 

M' 1 M! 2 ,M 1 M 2 



c 5 (M[M^,M 1 M 2 ;R)/R 5 
x \M[M 2 ^){M 1 M 2 \ (13) 



and 

h^ = h + J2 Ce ( M i M 2 ' MiMa ; ^l R& 

M[M^,MiM 2 

x |M{M^)(A/iAf 2 | (14) 

in the basis of the Zeeman substates \M\M<z) of the pair 
state nin2«. The components of the eigenvectors of the 
matrices C^(R) or Ce(R) are the coefficients of the ex- 
pansion of the eigenvectors of H^ or H ^ in this ba- 
sis. For the first order quadrupole interaction, these co- 
efficients are entirely determined by the angular factors 
D22 (o:M{ M 2 , olM\ M2; R) and are independent of ri\ and 
n%. They are given in table U for the cases of interest in 
this work. Because the C5 functions depend on S only 
through an overall factor, the eigenstates of H ( 5 > are the 
same for singlet and triplet states. They are identical to 
those obtained in Ref. [13[ for a spinless alkali atom. 

That f2 is a good quantum number when R — z origi- 
nates from the fact that Jtz commutes with W mt for this 
particular orientation of the interatomic axis. In contrast 
J t 2 does normally not commute with H' int . Nonetheless, 
as indicated in table HI some of the eigenstates of H^ 
for R = z are also eigenstates of J t 2 . In view of this 
fact, we label the simultaneous eigenstates of H^ and 
J t by a number K such that the corresponding eigen- 
values of J t 2 are K 2 K(K + 1). For eigenstates of H^ 
that are not eigenstates of J t 2 , we assign the number K 
to the eigenstate whose components in the Zeeman basis 
are closest to those of the eigenstate of J 2 with eigen- 
value Ti 2 K{K + 1) for the same fl. Doing so leads to the 
assignments indicated in table |U One may observe that 
the eigenstates are symmetric upon the interchange of 
the states of atom 1 and atom 2 for K even and antisym- 
metric for K odd. (The relevance of J 2 in this context 
had already been recognized in reference }12j|.) 

We label the eigenstates of H^ in the same fashion. 
As noted in table HI several of these eigenstates are also 
eigenstates of H^ and of J 2 . However, this is not the 
case in general. In particular, many of the eigenstates of 
ij( 6 ) vary with m and ri2, unlike the eigenstates of H^ 5 \ 
The correspondence between even and odd values of K 
and the symmetry under the interchange of the states of 
atoms 1 and 2 is nonetheless the same. 

In general, the composition of a pair state |AfiAf2) in 
terms of energy eigenstates will therefore depend both 
on the orientation of the internuclear axis and on R, al- 
though the dependence on R will be negligible when the 
expansion (J7J is completely dominated either by the 1/R 6 
term or by the 1/R 5 term. Due to the differences in the 
eigenvectors of H^ and of ff' 6 ', the total energy shift 
at the interatomic distances where \C§/R 5 \ ss \Cq/R 6 \ is 
best calculated by diagonalizing the Hamiltonian 

i? (5+6) = H Q + Yl \cs{M[M' 2 , M 1 M 2 ;R)/R 5 

M[M^MiM 2 

+c 6 (M[M^MiM 2 ;R)/R 6 } \M[M^}{MiM 2 \. (15) 
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TABLE I. The eigenstates of the _ff' 5 ' Hamiltonian in terms of the \M\, M 2 ) Zeeman substates of the pair states, for the case 
where the interatomic axis is aligned with the z-axis. It is assumed that L% — L2, Si — S2 and Ji = J2. The eigenstates listed 
in this table are the same as those given in Table 1 of Ref. [13( |. Note a: C5 = for this state. Note b: The state specified in the 
right-hand column is an eigenstate of J t and of H^ ' as well as of H^ ' for any n. Note c: The K — 0, fl — and K = 2, fl — 
eigenstates of H^ 5 ' are also eigenstates of J t 2 for J = 1; however, the K — 0, Q. = and K — 2, Q — eigenstates of H^ 6 ' are 
not. 



Obtaining the dispersion coefficients thus largely re- 
duces to a computation of radial matrix elements and of 
angular terms. We calculate the former using the an- 
alytical expressions derived in references [26H28I ] in the 
framework of the Coulomb approximation. This ap- 
proach yields accurate results for sufficiently high prin- 
cipal quantum numbers and does not require any other 
input than orbital angular momentum quantum numbers 
and binding energies. The latter are obtained from ex- 
perimental data. 

The calculation of the Cg coefficients also involves a 
summation over intermediate pair states. As is illus- 
trated by figure [T] the sum is dominated by the pair 
states with the smallest Forster defect A. Including in 
the sum the 15,000 (or thereabout) pair states with the 
smallest values of A was sufficient to ensure convergence 
of the Cq coefficients to four significant figures. (The 
intermediate states included were restricted to principal 
quantum numbers in the range 10 < n'{,ri2 < 100.) This 
rate of con verg ence is similar to that observed in the al- 
kali metals [13| . 

Before closing this section, we briefly comment on the 
applicability of the perturbative approach at the inter- 
atomic separations relevant for cold Rydberg gases ex- 
periments (typically 2 to 10 /mi). A comparison between 
a non-perturbative calculation of the energy shift due to 
the dipole-dipole interaction and the prediction of the 



leading-order perturbative calculation (AE = Cq/R 6 ) is 
shown in figure 2. The non-perturbative shift was calcu- 
lated by diagonalizing an Hamiltonian matrix of compo- 
nents 



-flpg — OpqOp 



Vpq, 



(16) 



where the indexes p and q run over all the pair states 
included in the calculation, £ p is the asymptotic energy 
of the pair state p, and 

V pq = Dn(a q Mi q , M 2q , a p Mi p M 2 p; R = z) 

x Rii{nipn2pap;ni q n 2q a q )/R 3 . (17) 

For the state considered in figure 2, the perturbative 
1/R 6 shift matches the non-perturbative result very well 
for R larger than about 1.5 /zm but there are large dif- 
ferences at smaller separations. (The barely visible dif- 
ferences noticeable at larger values of R originate from 
differences in the set of states taken into account: only 
4,000 pair states were included in the non-perturbative 
calculation.) Perturbation theory breaks down when AE 
is comparable to or exceeds the Forster defect with the 
nearest pair state, A, i.e., at a separation R np such that 
Ca/R^p ~ A. Hence R np scales with the principal quan- 
tum number approximately like n 7 ' 3 : increasing n from 
50 to 100 increases R np by about a factor 5, which may 
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FIG. 1. (colour online) Contribution of individual interme- 
diate 1 So pair states to the C& coefficient of the 5s50p 1 P\ 
state of strontium (top), and the inverse of the corresponding 
Forster defect (bottom). 
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FIG. 2. (colour online) Comparison between the energy shift 
calculated non-perturbatively (black dash-dotted curve), the 
Ce/R 6 energy shift (green solid curve) and the C5/R 5 en- 
ergy shift (red dashed curve) for a pair of strontium atoms, 
both in the 5s50p Pi state, at a distance R from each other. 
(The results shown are for the K — 2, Q — eigenstate, the 
only SI = pair state for which C5 is non-zero). The shaded 
area represents the region where R is smaller than the LeRoy 
radius. 



B. Energy Levels 

The calculation of the C5 and, particularly, the Cq 
coefficients requires an accurate knowledge of the bind- 
ing energies of all the relevant atomic states. However, 
the quality of the available experimental binding ener- 
gies varies from series to series, as some are known only 
through laser spectroscopy measurements while others 
have been measured using microwave spectroscopy. We 
have used the microwave measurements where available, 
in view of their normally higher degree of accuracy. The 
random errors on the energy levels obtained by laser spec- 
troscopy are often large enough to affect the Cq coeffi- 
cients significantly. To reduce their impact, we fit the 
corresponding energy levels to the Rydberg-Ritz formula 
for the quantum defect |2Jj. Thus, for a state of principal 
quantum number n, binding energy E\, and Rydberg con- 
stant i? a , we express E^ in terms of the quantum defect, 
6 = n - [i? a /£b] 1/2 , and write 



S = Sn 



preclude the use of the corresponding C5 and Cg coeffi- 
cients in systems where the typical interatomic spacing 
is a few microns. It is clear that nonperturbative ef- 
fects become increasingly important at high n and may 
be predominant at or near Forster resonances, where the 
Cq coefficient is large and the Forster defect A is small. 
(The energy shift due to the quadrupole-quadrupole in- 
teraction, represented by a dashed red curve in figure 2, 
will be discussed in Section ITTT1 ) 



(n-5 ) 2 (n-So) 4 



(18) 



We set 5k = for fc > 4 and obtain <5o, 82 and 84 by least- 
square fitting to the data. The resulting values of these 
coefficients are shown in tablelllltogether with the ranges 
of principal quantum numbers used in the fit. The un- 
certainties quoted in the tables were obtained by consid- 
ering the variation of the x 2 function about its minimum 
J30| . Correlations in these uncertainties, given by the off- 
diagonal components of the covariance matrix, are negli- 
gibly small. We used the experimental energies directly 



for strongly perturbed series that could not be fitted in 
this way and the published quantum defects for the series 
measured by microwave spectroscopy. 

In the case of strontium, no microwave measurements 
of energy levels are available but laser measurements have 
been made for a number of series over wide ranges of val- 
ues of n I22T . l3ll - [39| . We used the energies of references 
[22L 1311-34] in view of their higher accuracy. The values 
of the Sq, 82 and 84 coefficients are given in Table [II] 
The Rydberg-Ritz formula fits the data well for all the 
series, with the exceptions of the 3 P2 and 1 5o series for 
which substantial departures were found for high princi- 
pal quantum numbers. (In the case of the Sq states, the 
departure has been attributed to collisional shift with for- 
eign gas [3l|.) For these two series, the fit was restricted 
to the values of n over which the Rydberg-Ritz formula 
could match the data. The 1 £>2 and 3 I?2 series exhibit 
strong configuration mixing [22]; nevertheless the energy 
levels are well described by the Rydberg-Ritz formula for 
n > 20. 

Microwave spectroscopy measurements have provided 
very precise energy levels for calcium [40] , |41| and ytter- 
bium 42j . We have supplemented these results with laser 
spectrosopy measurements of the 3 P2 series of calcium 
[34j and of the 1 F 3 series of ytterbium [43|]. Altogether, 
though, fewer series have been measured for these ele- 
ments than for strontium (table HT)). Moreover, the 1 £>2 
series of calcium and the 1 F^ series of ytterbium, which 
are highly perturbed, cannot be fitted to the Rydberg- 
Ritz formula; as a consequence, we could not extrapo- 
lated the measured energies to other values of n. These 
limitations reduced the number of states for which we 
could calculate the Cq coefficient. 



Uncertainties 



Here, Z is the atomic number and B, a and /3 are pa- 
rameters obtained by least-squares fitting of the energy 
levels to the experimental data. The radial matrix el- 
ements were found to differ typically by about 0.1 a.u. 
between the two calculations, which translates to differ- 
ences in the values of the Cq coefficients of about 0.5% 
for 71 « 20 and less for larger values of n. These differ- 
ences are smaller than the uncertainty originating from 
the error in the energies for this series. However, it may 
be that the error introduced by the Coulomb approxi- 
mation dominates the total error where all the relevant 
energies are known accurately from microwave measure- 
ments. 



III. RESULTS AND DISCUSSION 

Tables of the C5 and Cq coefficients for the series 
listed in table IIIII are provided in the on-line supple- 
mentary data accompanying this paper and form the 
main results of this work. We only consider pair states 
where both atoms are in the same Rydberg state - 
thus n\ = n 2 = n, L\ = L 2 = L, Si = S2 = L and 
•h = h = J- Tables of the c 5 {M[M^, M 1 M 2 ; R) and 
the cq(M' 1 M! 2 , MiM 2 ; R) functions are also provided in 
the Supplementary Information for R = z, i.e., for the 
case where the internuclear axis is aligned with the axis 
of quantization of the angular momenta. These results 
can be used, e.g., to obtain the eigenvectors of H^ and 
H( 6 \ Simple polynomial fits to selected Cq coefficients 
are given in table IIVI 

An overview of these results is presented in the next 
sections. We first examine the form of the long-range 
interaction, then consider the C5 and Cq coefficients in 
more detail. 



Errors in the binding energies of the relevant states 
dominate the uncertainty on the dispersion coefficients 
for most series. The fit of the experimental quantum 
defects to the Rydberg-Ritz formula reduces this uncer- 
tainty in smoothing out the random scatter in the data. 
We estimated the uncertainty on the values of the C5 and 
Cq coefficients arising from the errors on the values of So , 
82 and 84 by varying each of these parameters one by one 
and adding the resulting differences in quadrature [30j j . 
(This procedure is likely to overestimate the total error 
for the series measured by microwave spectroscopy 40].) 
The error introduced by extrapolating the Rydberg-Ritz 
formula to outside the range of principal quantum num- 
bers used in the fits could not be ascertained. 

Another source of error is the use of the Coulomb ap- 
proximation to calculate the radial matrix elements. We 
have compared the Coulomb approximation for the 1 Sq 
series of strontium to a calculation of the Cq coefficients 
using a model potential of the form 

V model = ~[1 + (Z- l)e- ar + Bre-P r ] . (19) 



A. Long-range interactions 

The relative strength of the first-order quadrupole- 
quadrupole interaction (C5) and the second-order dipole- 
dipole interaction (Cq) is illustrated in figure[2J As is well 
known, the C5/R 5 term dominates the energy shift at 
large separation R (for the symmetries for which C5 7^ 0) , 
while the Cq/R 6 term is significant below a critical radius 
Re = Cq/C$. How large R c is depends on the symmetry 
of the eigenstate and on n, see figure [3] (R c is roughly 
proportional to n 3 ). For most states, the Cq/R 6 term 
dominates the energy shift up to interatomic distances 
at which this shift is too small to be relevant for ex- 
periments. However for certain symmetries, such as the 
fl, = 1 state represented in the figure, the second order 
dipole-dipole shift is unusually small due to a vanishing 
angular factor [12]. In this case, the quadrupole interac- 
tion may be of significant importance when considering 
a dipole blockade. 

Given that the quadrupole interaction can normally be 
neglected at the atomic densities involved in cold Ryd- 



Atom 


Series 


do 


5 2 


Si 


Fitted Range 


Ref. 


Sr 


5sns 1 So 


3.26896(2) 


-0.138(7) 


0.9 (6) 


14 < n < 34 


[31J 




5snp 'P\ 


2.7295(7) 


-4.67(4) 


-157 (2) 


10 < n < 29 


[33] 




5snd 'D 2 


2.3807(2) 


-39.41(6) 


-109 (2)xl0 1 


20 < n < 50 


[31] 




5snf 1 F 3 


0.089(1) 


-2.0(2) 


3 (2)xl0 1 


10 < n < 25 


[33] 




5sns 3 Si 


3.371(2) 


0.5(2) 


-1 (2)xl0 1 


13 < n < 45 


[32] 




5snp 3 P 2 


2.8719(2) 


0.446(5) 


-1.9 (1) 


8 < n < 18 


[34] 




5snp 3 Pi 


2.8824(2) 


0.407(5) 


-1.3 (1) 


8 < n < 22 


[34] 




bsnp Po 


2.8866(1) 


0.44(1) 


-1.9 (1) 


8 < n < 15 


[34] 




5snd 3 D 3 


2.63(1) 


-42.3(3) 


-18 (1) xlO 3 


20 < n < 45 


[32] 




hsnd 3 D 2 


2.636(5) 


-1(2) 


-9.8 (9)xl0 3 


22 < n < 37 


[22] 




5snd 3 Di 


2.658(6) 


3(2) 


-8.8 (7)xl0 3 


20 < n < 32 


[32] 




5snf S F 4 


0.120(1) 


-2.4(2) 


12 (2)xl0 1 


10 < n < 24 


[33] 




5snf 3 F 3 


0.120(1) 


-2.2(2) 


12 (2)xl0 1 


10 < n < 24 


[33] 




5snf 3 F 2 


0.120(1) 


-2.2(2) 


12 (2)xl0 1 


10 < n < 24 


[33J 


Ca 


4sns 'So 


2.337930 (3) 


-3.96 (10) 






[40] 




4snp 1 P 1 


1.885584 (3) 


-O.II4 (3) 


-23.8 (25) 




[40] 




Asnd 'D 2 




Highly perturbed 


series 




[40] 




4snf 'F 3 


0.09864 (9) 


-1.29 (9) 


36 




[41] 




Asns 3 Si 


2.440956 (3) 


0.350 (3) 






[40] 




Asnp 3 P 2 


1.9549(8) 


2.5(1) 


-16 (1) xlO 1 


12 < n < 60 


[34] 




Asnp 3 P 1 


1.964709 (3) 


0.228 (3) 






[40] 




Asnd 3 D 2 


0.8859(5) a 


0.13 (4) & 






[40] 




Asnd 3 Di 


0.8833(5) b 


-0.02 (4) b 






[40] 


Yb 


4/ 14 6sns ^o 


4.27914 (4) 


-7.06 (6) 


565 (25) 




[42] 




Af' 4 6snp 1 P 1 


3.95433 (5) 


-12.33 (6) 


1729 (27) 




[42] 




Af' 4 6snd'D 2 


2.71363 (4) 


-2.01 (4) 






[42] 




4/ 14 6m/ X F 3 




Highly perturbed 


series 




[43] 



TABLE II. The Rydberg-Ritz parameters for strontium, calcium and ytterbium used in the calculation of the C5 and Ce 
coefficients. Uncertainties in the last digits are given in brackets. The parameters printed in italic are quoted from the sources 
given in the last column of the table. The others were obtained in this work. For those, the last column gives the reference to the 
sources of the spectroscopic data used in the calculation and the fifth column the range of principal quantum numbers included 
in the fit. The Rydberg constants for strontium, calcium and ytterbium are Rs T = 109 736.627cm - [31[, Rc a — 109 735.81cm -1 
[U and 7? Y b = 109 736.96cm" 1 [3 respectively. 



a Due to a perturber, a term 9.08(9) x 10 
b Due to a perturber, a term 8.51(9) x 10~ 4 



[(n-S r 2 -0.01676700]" 
[(n -So)- 2 -0.01685410]' 



must be added to equation 
must be added to equation 



18) 
18) 



for this series, 
for this series. 



Atom 



Available C§ and ce coefficients 



Available C5 and C5 coefficients 



'So J 5i ^1 J P ,i,2 l D 2 ' 
'So ( J Pi) 3 Pi Cd 2 
'So 'Pi Cd 2 ) 



L Pl d Pl,2 l D 2 "ZJi.2,3 

'Pi 3 Pi, 2 Cd 2 ) 3 d 1 , 2 

'Pi 'D 2 



Strontium 

Calcium 

Ytterbium 



Di 



2.3 



TABLE III. Index of the coefficients tabulated in the supplementary data. Results are normally given for 30 < n < 70. 
However, only estimates for a reduced range of principal quantum numbers are given for the series indicated between brackets, 
due to a lack of spectroscopic data and the questionability of a single active electron treatment for this calculation. 



berg gases experiments, we consider the C5 coefficients 
only briefly in the following, before discussing the Cg co- 
efficients in greater detail. 



B. C5 coefficients 



As they depend on the energy levels only through 
radial matrix elements, the C5 coefficients exhibit less 
structure than the Cq coefficients, which also depend on 
the inverse of the Forster defect, 1/A. Due to selection 
rules on the orbital angular momentum, the first order 
quadrupole interaction vanishes for many of the series. 



Atom 



Series 



\M., 



5sns 1 So 
5sns 3 Si 
5snp x Pi 
5snd 1 D2 



a 

3.2xl0" d 
-2.387xl0~ 3 
-1.24xl0" 4 
-1.65xl0" 3 



Fractional Err or 



Sr 



4sns i 5'o 
4/ i4 6sns ^o 
4/ 14 6snp x Pi 



1.793x10^ 0.3190 

— — - — : — e : — _ _ _ — r 



0.51 


3.6 


1.211 


-21.18 


0.0349 


1.03 


0.365 


-7.05 



0.02 
0.01 
0.002 
0.02 



Ca 



-9.84xl0" 5 
-7.74xl0" 4 



-1.338 



0.0001 



Yb 



0.0234 
0.167 



-0.421 
-2.73 



0.003 
0.0009 



TABLE IV. Polynomial fits to the Ce coefficients for stretched states of strontium, calcium and ytterbium (Mi = M2 = 
Mj — ±J with respect to the internuclear axis). C& = n 11 (an 2 + bn + c) and the coefficients a, b and c are expressed in a.u. 
(Ce (GHz/im 6 ) = 1.4448 x 10~ 19 Ce (a.u.).) These fits are valid for 30 < n < 70. The fractional error quoted is the uncertainty 
in Cq for n — 50 due to the uncertainty in the energy levels. For all the series considered in the table this fractional error 
exceeds the deviation of the fitting polynomial from the calculated C§ coefficients. 
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FIG. 3. (colour online) The energy shift Cq/R 6 at the criti- 
cal radius R c where \C^/R 5 \ = |Ce/-R 6 |, for 5snp 1 Pi states 
of strontium. The principal quantum number, n, increases 
from 30 to 70 and from left to right in this figure. The total 
energy shift must be calculated by diagonalizing the Hamilto- 
nian denned by equation JTSJl where |C 5 /_R 5 | w \C e /R 6 \. The 
C5 coefficient is negative and the Cq coefficient is positive for 
Q — 1, which means that for this symmetry the total en- 
ergy shift vanishes at R w R c and the long-range interaction 
changes from repulsive below R c to attractive beyond R c . 



In particular, C5 = for J = 0. 

In general, the C5 coefficient scales with n like n 8 [13] ■ 
Apart for this scaling, these coefficients have a similar 
magnitude for most states. However, their sign varies 
from symmetry to symmetry. For the J = 1 states, the 
C5 coefficient is non-zero only for the K — 2 eigenstates. 
The corresponding energy shifts are generally one order 
of magnitude smaller for the K = 2, fi = ±2 eigenstates 
than for the K = 2, fl = ±1 and K = 2, Q = eigen- 
states, which have C5 coefficients closer in magnitude but 
opposite in sign. 

A large number of eigenstates have a non-vanishing 
quadrupole interaction for 3 — 1. The corresponding val- 
ues of C5 tend to arrange themselves evenly about zero, 



without marked differences between strontium, calcium 
and ytterbium. 

In general, most triplet states were found to have 
weaker quadrupole interactions than the singlet states 
for the same value of n, due to a smaller value of D22. 



C. Ca coefficients 

First of all, we consider the Rydberg series with or- 
bital angular momentum L = 0. In alkali metals, the S 
states are widely used in experiments as the interaction 
is repulsive, which reduces the effect of ionizing interac- 
tions, and is nearly independent of 9. In contrast to the 
alkali metals, there are two L = Rydberg series in two- 
electron atoms. The ^o states of the bosonic isotopes 
are particularly appealing for experiments as they have 
no magnetic sublevels. The Cq coefficients for the 1 Sq 
series of strontium, calcium, ytterbium and rubidium are 
compared in figure |H The three divalent atoms exhibit 
dramatically different behaviour. For Sr, the interaction 
is attractive]] while for Ca it is repulsive (although weaker 
than in Rb) . The interaction is also repulsive in the case 
of Yb, however the scaled Cq coefficients are over an or- 
der of magnitude smaller than for the other two species. 
These differences entirely arise from differences in the 
energy level spacing between these atoms. 

The isotropically attractive interaction for strontium 
may have uses in many-body entangled states [2(| and 
in non- linear self- focussing schemes [45[ . This interesting 
feature is not unique to this atom. We have estimated the 
Cq coefficient for the 1 Sq series of magnesium, mer cury 
and zinc using the quantum defects of references |46h49J | , 
and found the interaction to be repulsive for magnesium 
and mercury but attractive for zinc. 



1 Anomalous behaviour in the Sr 5sns 1 So energy levels 13111 for 
n > 36 (attributed to collisional shifts) leads to a discrepancy 
between the results in Table llVl and previous work |2(I . Exper- 
imental values were used in [201 ] . while we fit the energy levels 
at n < 34 where these effects are absent. The interaction is 
attractive in both cases. 
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FIG. 4. (colour online) The scaled C$ coefficient for the 2 5i/2 
series of rubidium [13|] and for the 1 Sq series of strontium, 
calcium and ytterbium. For the latter three atoms, the solid 
line represents the polynomial fit of Table IIVI The error bars 
are multiplied by 200 for calcium and ytterbium. 



Also of note is the weak interaction in the ytterbium 
So states, which illustrates the importance of performing 
detailed calculations for each series. Such a weak interac- 
tion would be detrimental to any experiment wishing to 
exploit the 1 So series of ytterbium to produce a Rydberg 
blockade. 

Whereas the dipole-dipole interaction is attractive in 
the 1 Sq series of strontium, it is repulsive in the 3 S\ series. 
For these states, H^ has six different eigenenergies for 
each n, three of which are doubly degenerate. (Different 
eigenstates have different energy shifts because of the fine 
structure of the intermediate 3 P states.) The resulting 
values of the scaled Cq coefficient, Cgn -11 , range from 
33.2 a.u. to 36.8 a.u. for n = 50. The strength of the in- 
teraction varies thus little between these eigenstates and 
is always larger than for the 2 S\/2 states of rubidium 
of same principal quantum number. These differences 
between eigenstates make the van der Waals interaction 
slightly anisotropic in the 3 S\ series, in that controlled 
excitation to a particular Mj state will generally excite 
a superposition of energy eigenstates whose composition 
will depend on the orientation of the internuclear axis of 
each pair of atoms in the cloud. Suppose, for example, 
that one prepares the two atoms of a pair in the same 
stretched state (Mi = M^ = 1). This state is an eigen- 
state of H^ only when the internuclear axis is aligned 
with the axis of quantization of the angular momenta. 
Otherwise, its composition in terms of eigenstates of H^ 
varies with the angle 9 between the two axes (figure O. 
However, as seen from the figure, the average Cq coeffi- 
cient varies by less than 10% between 8 = and 9 = it/2. 
This means that there are two nearly-isotropic S series 
in strontium with C% coefficients of opposite signs, and 
therefore that the interaction can be tailored to experi- 
mental requirements by modifying the excitation scheme 




FIG. 5. (colour online) Left: The expectation value of the 
scaled C§ coefficient for the n — 50 stretched states of the 3 Si 
(top) and 1 Pi (bottom) series of Sr, |1, 1), as a function of the 
angle 9 between the direction of quantization of the angular 
momenta and the internuclear axis. Right: The composition 
in energy eigenstates of these two stretched states. Red dash- 
dotted curves: K — 0, fi = 0. Green solid curves: K = 
2, Q — 0. Blue dotted curves: sum of K = 2, Q — 1 and 
K = 2, Q — — 1. Black dashed curves: sum of K — 2, £2 = 2 
and K = 2,9. = 2. The K = 1,0, = and K = 1,9. = 
±1 eigenstates are orthogonal to |1, 1) as they have opposite 
symmetry under the interchange of the states of atoms 1 and 
2. 



for the same atom. 

The difference in Cg coefficients between the different 
eigenstates of H^ is more significant in the 1 P\ series 
than in the 3 S\ series, leading to a larger relative varia- 
tion of the expectation value of the energy on the state 
| Mi = 1,M<2 = 1) (bottom row of figure O. The com- 
position of | Mi = 1, M2 = 1) in terms of the eigenstates 
of H ( 6 ' is very similar for the two series — it is actually 
identical for the £1 = ±1 and fi = ±2 eigenstates; how- 
ever, in the 3 S\ states it is governed by the orientation 
of the spin of the four valence electrons and in the X P\ 
states by the orientation of the two Rydberg p-orbitals. 

Figure [S] gives a snapshot of the strength and sign 
of the van der Waals interaction at n = 40. In the 
cases of the 1 £>2 and 3 D% series of Sr and of the l P\ 
series of Yb, the C@ coefficient is positive for some of 
the eigenstates and negative for the others, and for cer- 
tain values of K and £1 its sign depends on n. These 
sign changes occur through the Cg coefficient smoothly 
passing through zero as n varies. (The C§ coefficients 
are never exactly zero in our calculations, contrary to 
those discussed in Ref. fl2j, as we take into account all 
the intermediate angular channels.) The sign of Cg also 
changes with n in the 3 P\ and 3 D^ series of Sr, but in 
these two cases the change is abrupt and occurs at a 
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FIG. 6. (colour online) The scaled Ce coefficient at n = 40 
for all the series and energy eigenstates considered in this 
work. The red markers indicate the value of Ca in the 
stretched eigenstates (|fi| = 2 J). (Ce (GHz^m 6 ) = 1.4448 x 
10- 19 C 6 (a.u.).) 



Forster resonance, namely when the Forster defect of the 
dominant channel for the series changes sign and almost 
vanishes (figure [7]). The corresponding near degeneracies 
are between the 5s35p 3 Pi + 5s35p 3 Pi and 5s35s 3 Si + 
5s36s 3 Si pair states (A/2ir = 68 MHz) and between the 
5s37d 3 D 2 + 5s37d 3 D 2 and 5s34/ 3 F 3 + 5s35/ 3 F 3 pah- 
states (A/27T = 3 MHz), respectively. These two reso- 
nances also give rise to abnormally large values of the Cq 
coefficient in the 3 P\ and 3 D 2 series of Sr, as can be no- 
ticed in figure EJ We have not found Forster resonances 
in Ca or Yb, or in any other series of Sr for principal 
quantum numbers in the range 30 < n < 70. 
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FIG. 7. (colour online) The K = 2, |fi| = 2 scaled C 6 coeffi- 
cient for the Sr 5snp 3 Pi series (blue circles), and the Forster 
defect for the 2 x5snp 3 Pi — » 5sns 3 Si+5s(n+l)s 3 Si channel 
(red solid line). (C* 6 (GHz/im 6 ) = 1.4448 x KT 19 C 6 (a.u.).) 



IV. CONCLUSIONS 

We have presented perturbative calculations of the 
long-range interaction between calcium, strontium and 
ytterbium Rydberg atoms. Extensive tables of the C5 
(first order quadrupole-quadrupole) and Cq (second or- 
der dipole-dipole) coefficients are provided in the sup- 
plementary data. These calculations are based on the 
Coulomb approximation, and accurate binding energies 
are required as an input. The experimental energy level 
data currently available is summarized by the quantum 
defects listed in table QI] 

The 1/R 6 interaction is generally dominant in the 
range of interatomic separation important for experi- 
ments, although in channels where this interaction is un- 
usually weak the quadrupole-quadrupole 1/R 5 interac- 
tion may become significant. Comparing the Cg coeffi- 
cients for the same series revealed significant variations 
between the species, illustrating the fact that small differ- 
ences in energy level spacings may have dramatic effects 
on Rydberg interactions. For example, the 1/R 6 interac- 
tion for the l So series was found to be attractive for stron- 
tium, repulsive for calcium and nearly zero for ytterbium. 
The sign of the interaction can also change between dif- 
ferent spin symmetries, as, e.g., in the ^o and 3 Si series 
of strontium. Significant variation was also observed for 
different symmetries within the same atomic species, in 
particular in the presence of a Forster resonance. Such 
resonances are less common in two-electron atoms than 
in alkali metals; only two instances were found in this 
work, both in triplet states of strontium. 
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Since H[ nt does not couple singlet to triplet states, Si = 
S[ and £2 = S' 2 . Using standard summation rules and 
evalua ting the angular components of the matrix ele- 
ments [501 ] yields the familiar result 



Appendix A: Angular Factors 

Obtaining the C5 and C§ coefficients involves the 
calculation of angular matrix elements of the inter- 
action Hamiltonian H ln t between unperturbed states 
of the form given by equation ([l}. We collect 
the angular integrals together and define a coefficient 
D/ !l fe 2 (a'M{M2,aMiM2;i?) depending on the quantum 
numbers of the states, on the order of the multipole tran- 
sition considered, and on the orientation of the internu- 
clear axis. Here we present a derivation of these angular 
coefficients. 

Recall that we assume that the inner electron is 
in a s-orbital. Hence, l\ — mi { = 0, L = l , 
Ml = rni a and C hm L l m =1 and the angular part 

of the calculation reduces to that of the matrix elements 
(L'S'J'M'j\Y Kp \LSJMj). However, 

(L'S'J'Mj\Y k>p \LSJMj) = 

/ / / "m ei ,m a /0m. So ,m t , OMs,M' s 



m o > m i' Ms rn s' ,ms 



x (C£ 



SAI, 



J'M'j 



H s m So ) u Ui Slfs I>J S'M' S 



x {l' m' Vo \Yk, p \l Q mi ) 

EST^ r .lM., r rJ'M', 

/ ( ^ J l mi SMs^l' 



m[ SM S 



mi ,m l i Ms 

x (l' m' v \Yk,p\l mi ) . 



(Al) 



{L'S'J'M'j\Y Kp \LSJMj) 



( yio+i'o+i / (2^o + l)(2J+l)(2fc+l) 

r l' o r J'M'j f J k J' \ ( . 9 , 

x L/ i o 0fc0°JAfjfep 1 /' S I \ ' ^ ' 



Making use of these results and of equation ^ , we ex- 
press the matrix elements of the Hamiltonian H[ at as 
sums of radial terms multiplied by the coefficients 



D klk2 (a' M[M' 2 ,aM 1 M 2] R) 



'47r(2fci + 2fc 2 )!(2Lx + 1)(2Z, 2 + 1) 
(2/c 1 )!(2fc 2 )!(2fc 1 + 2fc 2 + l) 

,L',0 „LiO 



x V(2Ji + l)(2J 2 + l)C^ u fci0 C^ u fe20 

•h fci J{ 1 f J2 fc 2 J 2 
I^i o Li\ J 1 _L 2 o -L 2 



fci+fc 2 



k 2 

/ _ Z^ 2^ Y ki+k 2 ,p{R) 

p=-(k 1 +k 2 ) pi=— fcl P2 = -k 2 

JLML 



E 



„ /-<fei+fe2,P ^J'-iM', 



kipi,k 2 p 2 JiMi,fcipi J 2 M 2 ,k 2 p 2 



c 
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